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Abstract 

Kerr black holes with scalar hair are solutions of the Einstein-Klein-Gordon field equations describing 
regular (on and outside an event horizon), asymptotically flat black holes with scalar hair [l . These black 
holes interpolate continuously between the Kerr solution and rotating boson stars in D = 4 spacetime 
dimensions. Here we provide details on their construction, discussing properties of the ansatz, the 
field equations, the boundary conditions and the numerical strategy. Then, we present an overview 
of the parameter space of the solutions, and describe in detail the space-time structure of the black 
holes exterior geometry and of the scalar field for a sample of reference solutions. Phenomenological 
properties of potential astrophysical interest are also discussed, and the stability properties and possible 
generalizations are commented on. As supplementary material to this paper we make available numerical 
data files for the sample of reference solutions discussed, for public use. 
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1 Introduction 


There is strong observational evidence that extremely compact and massive objects populate the Universe. 
One piece of evidence comes from our own galactic centre, from the bright radio source Sagittarius A*. 
From the Keplerian orbits of stars in its vicinity, its mass has been estimated as 4.1 x 10 6 M© (where 
M©= 1 solar mass) [2] and a size constraint of 6 light hours has been placed on the object [3]. The best 
theoretical candidate from well established physical models which fits these observational data is a black 
hole (BH). Thus Sagittarius A*, as well as other similar compact objects at the center of other spiral and 
elliptic galaxies, is commonly referred to as a supermassive black hole. BH candidates of this type have been 
found within a mass range between 10 6 and 10 10 M© [I] and they are thought to play a central role in both 
the formation and growth of their host galaxies [5]; thus understanding them is of vital importance for the 
models of structure formation in the Universe. 

Another case of extremely compact and massive objects is found in binary systems in our galaxy, where 
strong X-ray sources exist. One example of such an object is Cygnus X-l, known since the 1960s 6 . The 
estimated masses of these 24 binaries were found to range between 5 and 30 M© [2]. For instance, Cygnus 
X-l has an estimated mass of 14.8 M© |7 . Neutron stars , as the most compact directly observable objects 
currently known, have masses which are below 3 M© [8] 9 . Thus, these 24 X-ray sources are thought to be 
BHs, which, within this mass range, are dubbed stellar mass black holes. 


At present, it is unknown if the BH candidates discussed in the previous two paragraphs are the paradig¬ 
matic BHs of general relativity, BHs as described by some alternative model, or even other types of compact 
objects without an event horizon. The next decade promises to shed light on both of these issues: we are on 
the verge of gathering observational evidence that will map the spacetime geometry close to these BH can¬ 
didates. This evidence will be obtained by both gravitational wave astronomy 
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and by large baseline 

interferometry measurements of the galactic center, using the Event Horizon Telescope (EHT). The latter, 
promises to resolve angular scales of the order of the horizon scale for the Sagittarius A* BH candidate [l3]. 


The EHT will study the so-called BH shadows 14 : the gravitational lensing and redshift effect due to the 


BH on the radiation from background sources, with respect to the observer. These forthcoming experiments 
make it particularly timely to explore alternative models to the general relativity BH paradigm, and their 
associated phenomenology 
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According to the general relativity BH paradigm the myriads of BHs that populate the Universe should, 
when near equilibrium, be well described by the Kerr metric p. The paradigm is supported on both a 


set of mathematical theorems - the uniqueness theorems 17 18] — but also on a conjecture - the no-hair 
conjecture [19 . The former established that, for vacuum Einstein’s equations, the only regular (on and 
outside an horizon) BH solution is Kerr. The latter extrapolates that, even for more generic forms of matter, 
the end-point of gravitational collapse should still be an exterior Kerr solution. Progress in this context was 
obtained for particular types of matter. One of the simplest types of “matter” often considered by physicists 
is provided by scalar fields. Since 2012, there is observational evidence that scalar fields exist in nature, 
by virtue of the discovery of a scalar particle at the Large Hadron Collider, at CERN, identified with the 
standard model Higgs boson [20,21 . But for decades, scalar fields have been considered in phenomenological 
models, in particular within gravitational physics. A notable example is cosmology, where various types of 
scalar fields have been used to model dark energy and dark matter. One reason is that scalar fields are well 
motivated by some high energy physics models such as string theory or scalar extensions of the standard 
model of particle physics. Yet another reason is that scalar fields may be considered as a proxy to realistic 
matter, such as some perfect fluids. 

It is therefore quite natural that in testing the no-hair idea, scalar fields were one of the first types 
of “matter” considered. This program was initiated by Chase 22 who established that “every zero-mass 


scalar field which is gravitationally coupled, static and asymptotically flat, becomes singular at a simply- 
connected event horizon”. In other words, a BH spacetime cannot support a regular massless scalar field in 
equilibrium with it; i.e. no BH (massless) scalar hair. Further “no-scalar-hair” theorems were developed by 
Bekenstein [23||24j who also considered massive scalar, vector and spin 2 fields (see also the review 25]) and 
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Hawking 26 , who showed that in the Brans-Dicke theory of gravity 27 , in which there is a scalar field non- 


minimally coupled to the geometry, the regular BH solutions are the same as in general relativity. Hawking’s 
theorem has actually been recently generalized [28] to more general scalar-tensor theories of gravity , of which 
the Brans-Dicke model is an early example^ 

A few remarks are in order concerning “no-scalar-hair” results for BHs. First of all, we are focusing 
on regular configurations satisfying the weak energy conditior0 Solutions where the scalar field diverges 
at the horizon are known (see e.g. [34]), but they appear to have no physical relevance. We are also 
considering independent scalar fields. Considering simultaneously gauge fields to which the scalar fields 
are non-minimally coupled leads to solutions 35,36 - notably p-brane type solutions in supergravity 37 


but these have no independent scalar charge and the scalar field vanishes when the electromagnetic field 
vanishes. Finally, we are focusing on four dimensional, asymptotically flat BHs. Considering, for instance, 
Anti-de-Sitter asymptotics can allow for hairy BHs, since the asymptotic nature of the spacetime yields a 
trapping mechanism from which the scalar field cannot escape. A thorough review of asymptotically flat 
BHs with scalar hair is given in 
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Considering more general kinds of matter, many counter-examples to the no-hair conjecture - at least 
in its weakest version, stating that no BH solutions with other forms of matter should exist, regardless 
of their stability - have actually been shown to exist. An influential pioneering example was built upon 
the Bartnik-Mckinnon “particle-like” solution of the Einstein-Yang-Mills equations 39 
it has been found that a BH could be added at the center of these solitons 
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Soon after that, 
yielding BHs with 

Yang-Mills hair also dubbed colored BHs. Other examples with a similar spirit, were obtained with other 
non-Abelian gauge fields, e.g. [42],43 ; see also the reviews 44p5 . Such counter-examples however, invariably 
use non-linear matter fields, and the resistance of the matter field against collapse into the BH is anchored to 
these non-linearities. Furthermore, many of these counter-examples use matter fields which, probably, have 
little astrophysical relevance, at macroscopic scales. As such, even though these examples show clearly the 
mathematical limitations of the no-hair idea, the question remains if there are astrophysically more realizable 
models of hairy BHs. 

Recently a new family of BHs with scalar hair was found [l] that is continuously connected to the Kerr 
family and yields a qualitatively new example of hairy BHs dubbed Kerr BHs with scalar hair (KBHsSH). 
Firstly because the scalar hair is not anchored on non-linear effects. The scalar field possesses a mass term 
but no self-interactions. Consequently, the hair can be seen in linearized theory, by considering the massive 
Klein-Gordon equation around the Kerr BH, and it was interpreted in [l] as a zero mode of the superradiant 
instability. We recall that the superradiant instability of a Kerr BH, in the presence of a massive scalar field, 
is a mechanism via which a field mode with frequency w and azimuthal harmonic index m is amplified when 
w < mQ#, where £2# is the Kerr BH horizon angular velocity 46 . By solving the Klein-Gordon equation 


on the Kerr background, real frequency bound states can be obtained when w = uiQh , corresponding 
to linearized (hence non-backreacting) hair, called scalar clouds [1,47- 50 . These are therefore, precisely 


between the decaying (w > mfln) and superradiant regimes, and thus can be faced as zero modes of the 
latter. Secondly, because they provide an example of scalar hair around rotating BHs, which, as discussed 
above has more potential to be astrophysically relevant. 

In this paper we will provide details of the construction and of the physical properties of this new type 
of hairy BHs, complementing the construction in |[l] and the discussions in 51 52 . The paper is organized 
as follows. In Section [2] we describe the Einstein-Klein-Gordon theory, the ansatz taken (complemented 
in Appendix 0 and exhibit the corresponding system of equations obtained. In Section [3] we discuss 
the boundary conditions to be imposed on the solutions (complemented in Appendix [Bj, introduce physical 
quantities and physical relations of interest, such as the first law of BH mechanics and Smarr relations and we 
discuss the numerical procedure that is used to obtain the solutions. In Section [4] we present the solutions, 
starting by describing relevant properties of the two limiting cases (boson stars (BSs) and a test scalar 


1 Scalar-tensor theo ries with higher derivatives (albeit with second order equations) can, however, accommodate spherically 
symmetric hairy BHs [29fj32 . 

2 Hair y BH solutions can be constructed by allowing a scalar potential which is not strictly positive, see e.g. the recent 
work |33l and the references therein. 
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field around Kerr BHs). In particular we shall present a set of five reference configurations that illustrate 
qualitatively distinct regions of the solution space and provide illustrative plots of these reference solutions (in 
Appendix [C| . In Section [ 5 ] we discuss physical properties of KBHsSH with astrophysical phenomenological 
interest, as the quadrupole moment and the orbital frequency at the innermost stable circular orbit. In 
Section [6] we address the issue of stability, exhibiting some of the known facts and the open questions. 
Finally, in Section [7] we enumerate a set of research directions and various generalizations of the solutions 
that can be addressed. 


2 Action, ansatz and equations of motion 

2.1 Action 

We shall be working with the Einstein-Klein-Gordon (EKG) field theory, describing a massive complex 
scalar field T minimally coupled to Einstein gravity. The model has the following action and corresponding 
EKG field equations, obtained from the variation of the action with respect to the metric and scalar field, 
respectively: 


5 = J d 4 xV=g ^^R- l -g ab {^ a ^, b + ^,a)-^*^ 


( 2 . 1 ) 


E ab = R ab - -g ab R - 8irG T ab = 0 , T ab = ’T a ’b b + 'I'*,/ 1 ', a - g ab 
□'b = /i 2 \b , 


-g cd {^ A + ** d tt, c ) + 


( 2 . 2 ) 


where G is Newton’s constant, /j the scalar field mass and T ab is the energy-momentum tensor of the scalar 
field. Observe that this model containing one complex scalar field is equivalent to a model with two real 
scalar fields. Indeed, writing 'b = 'b F + fib 7 , where \b F , 'b 1 are two real scalar fields, the action becomes 


S = j d 4 x^~g ^R - g ab (<b fi a <b F + - g 2 [(tf*) 2 + (^ / ) 2 ] 


(2.3) 


which describes Einstein’s gravity minimally coupled to two real massive scalar fields, with the same mass. 
The fact that one has two real scalar degrees of freedom is the reason why the solutions we describe circumvent 
the theorem in 


53 


The action (2.1) is invariant under the global 17(1) transformation T -A e m T, where a is constant. Thus, 


the scalar 4-current, j a = —i(T*9 a T — T9 a T*), is conserved: j a a = 0. It follows that integrating the 
timelike component of this 4-current in a spacelike slice E yields a conserved quantity - the Noether charge : 


Q -L f 


(2.4) 


At a microscopic level, this Noether charge counts the number of scalar particles. In the case an event 
horizon is present, Noether charge conservation does not prevent the scalar field from falling into the BH; 
rather, there is a continuity equation relating the decrease of the Noether charge with the scalar flux through 
the horizon. As such the scalar field may completely disappear through the horizon; moreover, since there 
is no Gauss law associated to the scalar field, it would leave no signature in the exterior spacetime. 


2.2 Ansatz 

Kerr BHs with scalar hair (KBHsSH) solutions are obtained by using the following ansatz for the metric and 
scalar field 


ds 2 = e 2Fl (^- + r 2 d6> 2 ) + e 2F2 r 2 sin 2 0(d<p - Wdt) 2 - e 2F °Ndt 2 , 


with N = 1-— , (2.5) 

r 
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( 2 . 6 ) 


where w is the scalar field frequency and m = ±1, ±2... is the azimuthal harmonic index; without loss of gen¬ 
erality, we take w > 0. The full configuration is therefore described by five functions of (r, 0): F 0 , Pi, F 2 , W, 4>. 
We observe that for m = 0 this is basically the ansatz used for obtaining rotating BSs [54|[55]. We further 
remark that, in the Kerr limit, this ansatz will yield the Kerr solution in a coordinate system which does not 
coincide with standard textbook coordinates for Kerr. The coordinate transformation from the coordinate 
system in (2.5) to Boyer-Lindquist coordinates is, however, quite simple. It is provided in Appendix [A[ 
Observe that KBHsSH are described by a metric ansatz with two Killing vector fields 


£ = 9 U 


and 


rj = d l 


<p m 


(2.7) 


£ and r] do not, however, generate symmetries of the full solution, since they do not preserve the expression 
of the scalar field. The only symmetry of the full solution is generated by the helicoidal vector field 


. w 

X = £+ —V: 
m 


( 2 . 8 ) 


since x ^ = 0- This combination is reminiscent of the null horizon generator for rotating BHs. KBHsSH will 
be obtained by choosing x to precisely coincide with such generator. 


Expanding a bit more on the two real fields picture, we remark that for T given by (2.6), it holds 
independently that x^ R = 0 and X^ 1 = 0, where 


= 0(r, 0) cos {rrup — wt ) , T 7 = 0(r, 0) sin (rrup — wt ) . 


(2.9) 


That the real and imaginary part of the complex scalar field are independently preserved by x implies that 
one real scalar field coinciding with either T R or T 7 can form stationary waves around a Kerr BH at linear 


level. These scalar clouds will be further discussed in Section [472] The existence of a fully non-linear solution 
with a stationary metric, however, requires the existence of two real scalar fields with opposite phases and 
the same mass, corresponding to one complex scalar field. The phase difference guarantee the cancellation 
of the t and p dependence in the total scalar energy-momentum tensor. 


2.3 Equations of motion 

Let us now address in detail the system of coupled partial differential equations (PDEs) obtained from this 


ansatz. Firstly, the explicit form of the Klein-Gordon (KG) equation (2.2) reads 


1 If rN' 

\rr + + 0,r(Eo ,r + ^2,r) + 0,0 + ^2,6>) + ( 1 + 


2N 


\ 2 cot 0 

) ^ 2 ' 10 ' ) 


e 2i?2 m 2 
r 2 sin 2 6 


c ^ F °(w — mW) 2 

-Sv- L + » 


N 


■ = 0 


This is a second order PDE for the function 0. Secondly, concerning the Einstein equations (2.2), the 
non-trivial components are E^ E*E%. These six equations are divided into two groups: four of 


these equations are solved together with the KG equation (2.10), yielding a coupled system of five PDEs on 


the five unknown functions. The remaining two Einstein equations are treated as constraints and used to 
check the numerical accuracy of the method. 


Each of the four Einstein equations we shall solve simultaneously with (2.10), should, as (2.10), have 


second derivatives of a single function. This is achieved by using the following combinations of the Einstein 
equations: 


K 

El 


Eh - El -El = 0 


■ E e e - E% 


El 


■ 2 WEI = 0 


( 2 . 11 ) 


El + E° g +E%- El - 2 WE* = 0 


% = 0 • 
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These four equations, multiplied by suitable factors, yield, respectively, second order equations for F 2 , Fp 
and W: 


1 ( 1 \ e 2F o+ 2 ^ 2 r 2 gin 2 0 / 2 

Fl ’" + ^N Fl ’ ee - \ F( ' r F 2,r + ^4,@4,d ^ ( *4 -T- 


1 11 _ tk _ Ttk 

r 2iV 


1 


rJV'N Fi, P cot6>F 0>(? / 2 

1 + W J 4 rVV + 8?rG ^ 4v' 


+ 


«2Fi 


TV 2 


e 2F °(w — rriW) 2 — 


e 2F2 m 2 N 
r 2 sin 2 0 


= 0 , 


11 1 e -2F 0 +2F 2r 2 gin 2 Q / 1 

F 2 ,rr + 4,00 + F 2 2 r + ^Fl e + 4,r4,r + ^4,04,0 + -^- (w 2 + W 2 0 


- Ft 


0 ,r 


cot0F o ,«\ , (, , rN'\ 3F 2 , r , 2cot9F 2 ,,, , 0 ^,e 2F ' ( 2 , 2e- 2F ’m 2 

-tf - ) + {■ 1+ w) — + -ar +8 '°-f (" + 7*f? 


= 0 , 


11 1 2F 0 +2F 2r 2 s i n 2 /3 / 1 

Fn.rr + -4 t4+ FL + ~FzF 2 a + F nr F 2r + -4r4 04 0 -—-- I 1C 2 + ~Fw 2 


r 2 N 


r 2 N 


r 2 N 


2N 


r 2 N 


3 rN'\ 2 4, r cot6>4 jfl 

+ 4iV 7 r + r 2 N 


+ 


N'F< 


2 ,r 


2N 


8nG- 


e 2Fl f2e~ 2F °{w-mWf 


N 


N 


4 r = o, 


w 


"T + ^2N W ’ de + ( 3i?2 - 

4 f TTr 3 cot QW 1 

rr + ^T 


4rN 

On the other hand, the two constraint 


,r ~ Fo, r )W :r + ~ Fq 

9Pi- 

32t rG 


equations are chosen 

4 -4 = 0, 


>,o)W, 

e 2F 1 -2F 2m ( w _ 
r 2 sin 2 6W 


-cj) 2 = 0 . 


to be 


and 


( 2 . 12 ) 


which yield, respectively, 


4 = 0, 


4 


,rr 


X 1 2 1 

^ F o,ee + F 2 ,rr - ^4,00 + 4,r 


- 2 ^Tl,r^2,r 
/ 3rN‘ 


2 N 


1 


- 72^4*4, 
Fp, r 1 


1 - 


^4% - 2 (4,r4,r - ^4,04,*) 

p —2F 0 +2F 2/r 2 s i n 2 /3 / 1 \ 1 

5 —k - 44+^ - 4^ 

. 2 cot # / 

ON ) &• 


rN'\ 


(2.13) 


- OF. 


= 0 , 


cmrl 


Fo,r0 + -^2,7*0 + ^O,r^O,0 + ^2,r^2,0 



— (^0,r^l,6> + ^l,r^0,6>) — (^l,r^2,6> + ^2,r^l,0) 

/4 n p—2F 0 +2F 2<r 2 s i n 2 Q 

- cot 0(Fi r - iV)-—-lT, r ^ 0 + 16FG0, r 0, 0 = 0 . 

r 27V 


In the next section we will address the boundary conditions and the numerical methods used to solve 
these equations. 
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3 Boundary conditions, quantities of interest and numerics 

KBHsSH are asymptotically flat solutions, which are regular on and outside an event horizon. In order to 
perform the numerical integration of the system of equations described in Section |2j appropriate boundary 
conditions must be imposed. In our study, we shall not consider the behaviour of the solutions inside the 
event horizon. The boundary conditions we have chosen implement asymptotic flatness and regularity at 
the horizon and at the symmetry axis. Let us describe these boundary conditions in detail. 

3.1 Boundary conditions 

Asymptotic boundary conditions. For the solutions to approach, at spatial infinity described by 
r —>> oo, a Minkowski spacetime background we require 


lim Fi = lim W = lim 0 = 0. (3.1) 

r—> oo r—toc r —>-00 

For any input parameters, one can obtain an asymptotic expression of the solution, compatible with these 
boundary conditions. This expression is given in given in Appendix |B} 

Axis boundary conditions. Axial symmetry and regularity impose 
on the symmetry axis, i.e. at 6 = 0, i r: 

doFi = d e W = 0 = 0. 

Moreover, the absence of conical singularities implies also that 

F!=F 2 , (3.3) 


the following boundary conditions 

(3.2) 


on the symmetry axis. 

Also, all solutions discussed in this work are symmetric w.r.t. a reflection on the equatorial plan^] As a 
result, it is enough to consider the range 0 < 0 < 7t/2 for the angular variable, the functions F0 W and 0; 
these satisfy the following boundary conditions on the equatorial plane 

d » F <\e =*/2 = d » W \e=.l 2 = M*=,r /2 = °- ( 3 - 4 ) 


Event horizon boundary conditions. The event horizon is located at a surface with constant radial 
variable r = r# > 0. The boundary conditions there and also the numerical treatment of the problem are 
simplified by introducing a new radial coordinate 


*=* ~r% 

Then a power series expansion near the horizon yields 


Fi = F^\6) + x 2 f[ 2 \o) + 0(x 4 ), 
W = tt H + 0(x 2 ), 


(3.5) 


(3.6) 

(3.7) 


and 


cfr = M0) + O(x 2 ), 


(3.8) 


where the constant Qh > 0 is shown to be the horizon angular velocity, see (3.14). The field equations 
together with (3.6)-(3.8) imply that this quantity obeys the condition 


w = rriVtu- (3.9) 

3 We have also found solutions with an anti-symmetric scalar field w.r.t. reflections along the equatorial plane, while the 

metric functions are still even parity. Such configurations, however, are hard to study systematically and are likely to be more 
unstable. 








This guarantees that the null geodesic generators of the horizon are tangent to the Killing vector field x, 
defined in (2.8). The physical significance of such identification is that there is no flux of scalar field into 
the BH, 

X^T = 0, (3.10) 

which is central to the existence of regular BHs with a stationary geometry and scalar hair. As discussed 
in the Introduction, the condition (3.9) is also related to the superradiance phenomenon. Furthermore, we 
note that the Einstein equation E e r = 0 implies that the difference F\ — Fq is constant on the horizon. In 
our scheme, however, we do not impose this condition, but rather use it as another test of the numerical 
accuracy of the solutions. 

To summarize, the boundary conditions at the horizon are 


d x FA = d x (j)\ = 0, 

d, L \ r — rH x r | r — r . 


w\ 


w 

m 


(3.11) 


3.2 Quantities of interest and a Smarr relation 

Most of the quantities of interest are encoded in the expression for the metric functions at the horizon or at 
infinity. Considering first horizon quantities, we introduce the Hawking temperature Th = ft/(27r), where 
k is the surface gravity defined as k 2 = — aXb)iy a ^X b )\r H > an d the event horizon area Ah of KBHsSH. 
These are computed as 


T h = 


Attth 


„F™(e)-F< 2 \e) 


Ah = 27r r 


H 


f 


dO g 


l0 e F[ 2 \e)+F?\e) 


(3.12) 

(3.13) 


Moreover, S = An j(iG) is, as usual, the BH entropy. Also, the event horizon velocity tin is fixed by the 
horizon value of the metric function W, 


o _ 

— ~ t- —- 


Qyt 

9tt 


= W 


(3.14) 


The ADM mass M and the angular momentum J are read from the asymptotic sub-leading behaviour 
of the metric functions: 


g tt = —e 2F °N + e 2F2 W 2 r 2 sin 2 9 = -1 + 


2 GM 


+ ..., g vt = — e 2F2 Wr 2 sin 2 9 = — 


2GJ 


1 2 9 


(3.15) 


As usual in (asymptotically flat) BH mechanics, the temperature, entropy and the global charges are 
related through a Smarr mass formula [56 , which for the KBHsSH reads 


M = 2 T h S + 2 n H (J - mQ) + M*, 


(3.16) 


where M*, given by 

/* nOO [*7 T 

-M^= dS a (2T^ b -T£ a )=4ir dr 

J S J rn ^ 0 

is the scalar field energy outside the BH, and 

- f dr f 

J r h J 0 


2 a a F 0 +2F 1+ F 2 (p 2 _ 2e -2 F2 w(w-mlK) ^ 


dO r 2 sin 0 


N 


OO r-7T 

Q = 47 t j dr I dd r 2 sin0 e~ F ° 
r H 


-F 0 +2Fi+F 2 m ( w - mW) ,2 

N V 


4> 2 (3.17) 


(3.18) 


is the conserved Noether charge (2.4). 
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Some of these physical quantities are also connected via the first law 

dM = TtfdS H- QjjjdJ. (3.19) 

A natural question in the context of the hairy BHs we are discussing is how much ADM energy and 
angular momentum is in the BH and how much is in the scalar hair outside the BH. This question can be 
addressed by noting that the ADM quantities M and J can be expressed as 

M M^ T Mh, J = mQ T Jhi (3.20) 


where Mh and Jh are the horizon mass and angular momentum, computed as Komar integrals, which, 
from (3.16) and (3.20), satisfy the relation 


Mh = 2 ThS + 2QhJr- 


(3.21) 


3.3 Numerical implementation 


As usual when dealing with gravitating massive scalar fields, the numerical integration is performed with 
dimensionless variables introduced by using natural units set by fi and G, 


r —> r/i, (/)—> (j)Mpi/y/ Tzr, w—tw/fi, 


(3.22) 


where Mj> z = G~ l is the Planck mass. As a result, the dependence on both G and /i disappears from the 
equations. Also, the global charges and all other quantities of interest are expressed in units set by g and G 
(note that, in order to simplify the output, we set G = 1 in what follows). 

In our approach, the EKG equations reduce to a set of five coupled non-linear elliptic partial differential 


equations for the functions T = (To, T\, T 2 , W\ 0), which are displayed in Section 2.3 These equations have 
been solved numerically subject to the boundary conditions introduced above. An important i ssue here 
concerns the status of the two constraint equations E r e = 0, E r r — Eq = 0, also presented in Section 


is addressed following Ref. 


57 


2.3 


This 

One notes that the Bianchi identities V M TV r = 0 and V ^E ^ 9 = 0, imply 
the Cauchy-Riemann relations dfEz + deTZi = 0, dfEi — 80^2 = 0, where we have defined 
IZ 2 = y/Mjry/N ( E r r — E ®)/2 and the new variable dr = Therefore the weighted constraints still satisfy 

Laplace-type equations. Then they are obeyed when one of them is satisfied on the boundary and the other 
at a single point 57 . From the boundary conditions we are imposing, it turns out that this is indeed the 


case, i.e. the numerical scheme is consistent. 

Our numerical treatment can be summarized as follows. The first step is to introduce a new radial variable 
x = x/(l + x) which maps the semi-infinite region [0, 00 ) to the finite region [0,1] [we recall x = yG 2 — 
with r the radial variable in the line element (2.5)]. This involves the following substitutions in the differential 
equations 


xTx —(1 - x) x 2 F, xx —> (1 - x) 2 T xx - 2(1 - x)7^ x . (3.23) 

Next, the equations for T are discretized on a grid in x and 0. Various grid choices have been considered, 
the number of grid points ranging between 280 x 20 and 90 x 70. Most of the results in this work have 
actually been found for an equidistant grid with 250 x 30 points. The grid covers the integration region 
0 < x < 1 and 0 < 0 < tt/2 . 

All numerical calculations have been performed by using the professional package fidisol/cadsol [58] , 
which uses a finite difference method with an arbitrary grid and arbitrary consistency order. This package 
has been extensively tested in the past by recovering numerous exact solutions in general relativity and 
field theory. Furthermore, some of the new solutions derived by using fidisol/cadsol were rederived 
subsequently by other groups with different numerical methods. 

This code requests the system of nonlinear partial differential equations to be written in the form 
P(x, 0 ; Tx, Fo', FxO, Fxx, Foe) = 0, subject to a set of boundary conditions on a rectangular domain. 
Besides that, fidisol/cadsol requires the Jacobian matrices for the equations w.r.t. the functions T and 
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their first and second derivatives, the boundary conditions, as well as some initial guess for the functions T. 
Indeed, this solver uses a Newton-Raphson method, which requires a good first guess in order to start a suc¬ 
cessful iteration procedure. Also, this software package provides error estimates for each unknown function, 
which allows judging the quality of the computed solution. The numerical error for the solutions reported in 
this work is estimated to be typically < 10 -3 . However, errors increase dramatically when studying solutions 
close to central inspiralling region of the BSs curve, see Figure |4| below (left panel). A detailed description 


of the numerical method and explicit examples are provided in 58 


As a further check of numerics, we have verified that the families of solutions with a varying frequency 
satisfy with a very good accuracy the first law of thermodynamics ( 3.19| ) and also the Smarr relation (3.16). 

In the scheme we have used, there are three input parameters: i) the frequency w and ii) the winding 
number m in the ansatz (2.6) for the scalar field T, together with iii) the event horizon radius r# in the 
metric ansatz ( |2.5| ). The number of nodes n of \<fi\ on the equatorial plane, as well as all other quantities 
of interest ( e.g . mass, angular momentum, Noether charge, Hawking temperature and horizon area) are 
computed from the numerical solution. Both here and in [l], for simplicity, we have restricted our study 
to fundamental configurations, i.e. with a nodeless scalar field, n = 0. Excited states are more difficult to 
investigate systematically and, in any case, are expected to be more unstable. Also, we have studied in a 
systematic way the BH solutions with m = 1; sets of solutions with m = 2,3,4 have also been constructed 
for several fixed values of w. 

In some of the calculations, we interpolate the resulting configurations on points between the chosen grid 
points, and then use these for a new guess on a finer grid. Finally, the compilation of the numerical output 
is done by using the software mathematic A. 

Let us close this section with some technical details on the systematic procedure we have used to scan 
the parameter space of the KBHsSH solutions. We start with a spinning BS solution - i.e. having r# = 0 - 
with given m, re, as initial guess for a KBHSH with a small event horizon radius. Then we slowly increase 
the value of r#, keeping m,w fixed (thus also the event horizon velocity f£#). The iterations converge, and, 
in principle, repeating the procedure we obtain in this way solutions for higher and higher values of r#. 
Then, for the coordinate system we are using, a maximal value of r# is a approached and a second branch of 
solutions emerges, extending backwards in r#. The limiting behaviour of this secondary branch of KBHsSH 
depends on the value of the input parameters re, m. We find numerical evidence for the existence of three 
different possible limiting configurations on this secondary branch; they can be: i) another BS solution with 
Tfj = 0; ii) an extremal KBHSH (which is also approached as r# -A 0), and, iii) a special set of vacuum 
Kerr solutions with r# > 0. Thus, in principle, the full set of KBHsSH with a given m can be constructed 
in this way, by repeating this procedure for different values of the scalar field frequency w. The results 
reported in this work are obtained from around three thousand solution points. For all these solutions we 
have monitored the Ricci and the Kretschmann scalars, and, at the level of the numerical accuracy, we have 
not observed any sign of a singular behaviour. 


4 KBHsSH: exterior space-time and scalar field structure 

The solutions obtained by solving the equations exhibited in Section [2] with the boundary conditions and 
numerical method discussed in Section [3] form a 5-parameter family. Three of these parameters are continu¬ 
ous: i) the ADM mass M, ii) the ADM angular momentum J and iii) the Noether charge Q. The remaining 
two parameters are discrete: iv) the azimuthal hamonic index m G Z, appearing in (2.6) and v) the number 
of nodes (along the equatorial plane) of the scalar function 0(r, 0), n G No. Observe that the latter does 
not appear explicitly in the ansatz. Observe also that this 5-parameter family contains both horizonless 
solutions - BSs - and hairless solutions - Kerr BHs. To describe these two limiting solutions it is convenient 
to introduce the normalized Noether charge 

9 = ^j' ^ 

Then, the domain of solutions consists of q G [0,1], with Kerr BHs corresponding to q = 0 and BSs to q = 1. 
This latter statement will be expanded below. To contextualize KBHsSH let us start by introducing two 
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ingredients: i) the BSs that arise in the q = 1 limit of this family; and ii) a test scalar field in the Kerr 
geometry, which arises in the q = 0 limit. 


4.1 Solitonic limit: boson stars (q — 1) 


In this case r# = 0 and the horizon is replaced with a regular origin. The corresponding solutions are well 
known in the literature - BSs arguably the physically most interesting gravitating solitons. BSs have 
found a variety of applications, being considered as possible BH mimickers and dark matter candidates - see 
e.g., for a recent review 59 . 

The study of BSs was initiated by the work of Kaup 60 and Rufhni and Bonazzala 61 more than 40 
years ago. They found globally regular, asymptotically flat, equilibrium solutions of the Einstein equations 
coupled with a massive complex scalar field, providing an explicit realization of Wheeler’s geons 62 in EKG 


theory. The BSs were regarded as ‘macroscopic quantum states\ which are prevented from gravitationally 
collapsing by Heisenberg’s uncertainty principle. 

The BSs studied in the pioneering works [60,61 are spherically symmetric. Rotating generalizations, 
were first studied in the work of Schunck and Mielke 54 where such configurations were constructed in the 
weakly relativistic regime for a large range of winding numbers m. Highly relativistic spinning BSs have 
been constructed for the first time by Yoshida and Eriguchi 55 for winding numbers m = 1,2. These 


results have been extended recently 63 , in particular by constructing higher winding number solutions. For 


completeness, let us mention that spinning BSs with a self-interacting potential of the Q-ball type (see the 
discussion of Section [7]) have been studied in 64,65 , where odd-parity solutions were also first addressed. 


In the present context, rotating BSs are a particular limit of KBHsSH; as such they form part of the 
boundary of their domain of existence. Therefore we have performed an independent study of their properties 
for m = 1,2,3, by using the same methods described above, i.e. by solving a boundary value problem for 


J 7 , with two input parameters: w and m. The ansatz is again (2.5) (with r# = 0) together with (2.6). The 
boundary conditions at the origin, r = 0 reacQ 


d r F 1 = 8 r F 2 = d r F 0 = d r W = 0 = 0. 


(4.2) 


The boundary conditions as r —)> oo and at 9 = 0, tt are similar to those that apply to the KBHsSH case and 
that were described in Section [3l 

In Figure [T| we illustrate the BSs solutions. In the left panel we can see the ADM mass M vs. the 




Figure 1: Rotating BSs with m = 1,2, 3 and n = 0 in an ADM mass vs. scalar field frequency diagram 
(left panel) and an ADM mass versus ADM angular momentum diagram (right panel). 

frequency w distribution for BSs with m = 1,2,3, represented by the solid red, dashed green and dotted 

4 Note that, however, most of the BS solutions have been computed with W = W/r. 
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blue lines, respectively. We focus on nodeless solutions, n = 0, since these are typical the most stable ones, 
and we use natural units set by fi. On this plot we see that BSs exist for re < fi\ this is a bound state 
condition. As we decrease the frequency the mass increases until a maximum value. This value is of the 
order of 1 //i (or Mj> z //i, reintroducing the Planck mass Mpi). Thus in order to have such a BS with the 


eV. While such 
-, more reasonable masses 


mass of the Sun we would need extremely light scalar particles, with masses around 10 -11 
light scalars have been suggested in string compactifications - the Axiverse 

lead to so called mini-BSs [67]. Further decreasing w one finds a minimal frequency (seen for both m = 1 
and 2; for m = 3 the data collected did not reach the minimal value of w but we expect a similar pattern) 
below which no BS solutions are found. The BS curve then seems to spiral towards a central region of the 
diagram where numerics become increasingly challenging. Qualitatively, this is also the behaviour found for 
spherically symmetric BSs, m — 0, in which case, a detailed investigation of the inspiraling behaviour was 
possible. 

In the right panel of Figure [l] we exhibit the ADM mass M vs. the ADM angular momentum J distri¬ 
bution for the same BSs. This curve zig-zags, as can be seen for m = 1 and partly for m = 2. Each branch 
of this zig-zag pattern corresponds to a branch of the curve on the left panel where the mass increases or 
decreases. In Section |T3| we will see how KBHsSH fit in these two diagrams. 

Let us remark on the conserved Noether charge carried by the BSs. As observed above, this cannot be 
transformed into a flux at infinity; it is given by an integral over a space-like slice of the time component 
of the 4-current. A simple calculation shows that the angular momentum carried by rotating BSs relates to 
the Noether charge as J = mQ 54,55 64 . Thus q = 1 for BSs as advertised above. This will be useful in 
parameterizing KBHsSH. 

Finally, one may wonder how ‘compact’ these BSs are. Since BSs have no surface, i.e. the scalar field 
decays exponentially towards infinity, cf. (B.l), there is no unique definition of the BS’s ‘radius’. One 
estimate is provided in the following way. Firstly, we note that the ‘perimeteral’ radius, i.e., a radial 
coordinate R such that a cir cumference along the equatorial plane has perimeter ~ 2irR, is related to the 
radial coordinate used in (2.5) as R = e F2 r. Secondly, we compute Rgg, the perimeteral radius containing 99% 
of the BS mass, Mgg. Finally, we define the inverse compactness by comparing Rgg with the Schwarzschild 
radius associated to 99% of the BS’s mass, Rschw = 2Mgg 


68 


Compactness = 


R 


'99 


2Mgg 


(4.3) 


The result for the inverse compactness of BSs with m m 1 is exhibited in Figure [2] With this measure, the 
inverse compactness is always greater than unity; in other words, BSs are less compact than BHs, as one 
would expect. 

As can be seen from Figure [2] the least compact BSs are the ones closer to vacuum, which can be orders 
of magnitude less compact than a BH with the same total mass. At maximal mass, BSs are already only a 
few times less compact than a BH and during the whole spiral remain very close to the compactness of a 
BH. In the inset of the figure one can also see the maximal value of the scalar field along the BS line. For 
each solution this value occurs at a different radial coordinate. Observe that </> max increases monotonically 
along the BS curve; thus it is non-degenerate and could be used to label uniquely BS solutions. Indeed, 
historically, spherically symmetric BS solutions were labelled by the value of 0 at the origin. For rotating 
BS solutions this value is zero, but </> max could be used as a label. 


4.2 Kerr limit: scalar clouds (q = 0) 


In this case T is ‘small’ and we can linearize the coupled EKG equations (2.2) in T; these become simply 
R^ v = 0 and DT = /i 2 T. Thus we are led to study the KG equation for a complex massive scalar field on 
a vacuum solution of the Einstein equations. We consider the Kerr solution. As we shall see, a particular 
subset of Kerr solutions defines another limit of KBHsSH and form another part of the boundary of their 
domain of existence. 

Consider the KG equation on the Kerr background in Boyer-Lindquist coordinates and separate variables 
in the standard fashion: T = Si m (6)f (r), where Sim are spheroidal harmonics. One obtains a 
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Figure 2: Inverse compactness of BSs with mm 1. The inset shows the maximal value of the scalar field 
along the BS line. 


linear, second order, ordinary differential equation for /(r) 1,69-72 . One may then search for bound-state 
type solutions of this equation. Requiring that the radial function decays exponentially with r, however, 
yields generically quasi-bound states which have a complex, rather than real, frequency. The imaginary part 
of the frequency manifests that, in general, the scalar field cannot be in equilibrium with the BH; in fact one 
would expect it to fall into the BH, and this is precisely what one finds if the background is a Schwarzschild 
BH. In this case the imaginary part of the frequency is always negative corresponding to the scalar field 
decaying into the BH. In the Kerr case, however, it turns out that there is a critical frequency given by the 
product of the azimuthal harmonic index, m, and the horizon angular velocity, fin'- Wc — rnfln- It defines 3 
qualitatively different cases. If the real part of the frequency is larger than the critical frequency, 7 Z(w) > w c , 
then the quasi-bound state decays with time. This is the behaviour described above for Schwarzschild and it 
is the typical behavior expected due to the purely ingoing boundary condition at the horizon. If 7 Z(w) < w c , 
however, the quasi-bound state grows in time, signaling an instability. This is the superradiant instability 



In this case the existence of scalar clouds yields a quantization condition determining one physical possible 
(physical) value of the BH mass, which determines the angular momentum. For the non-extremal case, the 
clouds with a fixed (n, Z,m) exist for a 1-parameter subspace of the 2-dimensional Kerr parameter space. 
This is shown in Figure [3j where we exhibit a mass vs. horizon angular velocity diagram for Kerr BHs. The 
latter exist below the black line which corresponds to extremal Kerr. The blue dotted lines correspond to 
the backgrounds that support scalar clouds with n = 0 and m = /, for m = 1, 2, 3,4 and 10. 

It is important to observe that, fixing m and n, the scalar cloud whose existence line occurs for smaller 
angular velocity - for fixed M - is the line with m = l [l 50 . This line divides the Kerr parameter space 
into Kerr BHs which are unstable against some scalar modes with that value of m and Kerr BHs that are 
stable against all scalar modes with that value of m. Fixing m, therefore, this line, plays a special role and 
we shall call it the fundamental existence line for the azimuthal harmonic index m. 
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Figure 3: Existence lines for scalar clouds around Kerr BHs. The inset shows the radial profile along the 
equator for two clouds with n m 0, m = l «= 1. The colors correspond to the points with the same color in 
the existence line. The clouds are regular at the horizon, they attain a maximum and decay exponentially 
towards infinity. From [I]. 


4.3 KBHsSH solutions (0 < q < 1) 

4.3.1 Overview 


Let us start with an overview of the full space of KBHsSH solutions. In Figure [4] we fill in the plots shown in 
Figure [l] with the domain of existence of KBHsSH - the shaded blue region, for solutions with n = 0, m = 1. 
Let us mention that here and in figures below, the domain was obtained by extrapolating to the continuum 
the results from discrete sets of (thousands of) numerical solutions, cf. Section 3.3 This can safely be done 
for most of the parameter space; however, we do not exclude a more complicated picture for a small region 
around the center of the BS spiral, which is more difficult to explore numerically. 




Figure 4: Domain of existence for KBHsSH with m = 1 and n = 0 (shaded blue region) in an ADM mass vs. 
scalar field frequency diagram (left panel) and an ADM mass vs. ADM angular momentum diagram (right 
panel). Improved from l]; note, in particular that the extremal KBHsSH line contains more data then in [l]. 

In the left panel of Figure [4] we see the domain of existence in the ADM mass vs. scalar field frequency 
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diagram. The region where KBHsSH exist is delimited by: 

i) the BS curve already discussed in Figure [l]- where q = 1; 

ii) by the subset of Kerr solutions that support the fundamental existence line of scalar clouds with 
n = 0, m = 1 - where q = 0. In particular, this demonstrates that these hairy BHs are the non-linear 
realization of scalar clouds; 

iii) and by an yet unseen (green dashed) curve, corresponding to extremal (i.e. zero temperature) KBHsSH. 
We have evidence, but not a complete analysis that along this line q varies from zero to one and that it 
is the end point of all constant q lines, which inspiral in a similar way to the BS line. We shall further 
comment on these extremal solutions in Section 14.3.41 

Based on the existing numerical data, we are confident that the same pattern for the domain of existence 
of KBHsSH occurs for other values of m. In Figure [l] (see also [l ) the BS lines for m = 1,2,3 are also 
exhibited. 

In the right panel of Figure [4] we return to the ADM mass vs. angular momentum plot. In Figure [l] only 
the red solid curve had been shown, corresponding to BSs with m = 1. Kerr BHs exist in the upper part 
of the diagram, above the black solid line, which corresponds to extremal Kerr and KBHsSH exist in the 
blue shaded area. As before, the dotted blue line is the Kerr limit, corresponding to q = 0. Three general 
observations can be made: 


r) The first observation is that KBHsSH can violate the Kerr bound, since there are solutions below 
the black solid line. This is not surprising, since it is known that BSs can violate this bound 78 


Since KBHsSH are continuously connected to BSs one would expect the same to occur - at least some 
solutions with q close to one - and that is precisely what we see in the plot. But one may wonder if 
the Kerr bound is still violated in terms of the horizon angular velocity and mass. This question turns 


out to have interesting implications and will be discussed in detail elsewhere 79 


b) The second observation is that there are KBHsSH with the same mass and angular momentum as 
Kerr BHs. In this sense, and because M, J are the only asymptotic charges, there is non-uniqueness. 
Further specifying q , however, seems to completely raise the degeneracy. At least we found no evidence 
that there are two distinct solutions with the same (M, J, q). 

c) The third observation is that, in contrast to the vacuum case, KBHsSH do not possess a static limit Q 
A lower bound for the horizon angular velocity is set by the minimal value of the scalar field frequency 
of the corresponding BSs. Also, the maximal value of the BSs mass and angular momentum set an 
upper bound for the global charges of KBHsSH. 


In the region of non-uniqueness one can compare the area or entropy of the Kerr BHs and KBHsSH with 
the same M, J [I]. One observes that in the common region KBHsSH are entropically favoured which means 
that they cannot decay adiabatically to Kerr BHs - Figure [5] (left panel) . 

In Figure [6] (left panel) we exhibit again the ADM mass vs. scalar field diagram, but now showing 
constant Q lines. We recall that Q is a measure of the amount of scalar field outside the BH. Comparing 
with Figure [4] (left panel) we observe that constant Q lines span the domain of existence of KBHsSH in a 
very different way to constant q lines, which inspiral in a similar way to the BS curve; the former interpolate 
between two BS solutions and roughly, more massive KBHsSH have also larger Q. Note that here the mass 
is the ADM mass, as such taking into account both the BH (horizon) mass and the energy of the scalar 
field outside the horizon, cf. (|3.20|). In Figure [6] (right panel) we exhibit lines of constant q KBHsSH in a 
J/M 2 vs. ADM mass diagram. In this ‘phase space’ KBHsSH can coexist with Kerr BHs (which exist in the 
pale blue shaded region) on either side of the fundamental existence line (q = 0), but lower (larger) mass 
solutions than the Kerr limiting solution occur for small (large) q. 


5 Which is in agreement with the theorem in 80 . 
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Figure 5: (Left panel) Domain of existence of KBHsSH (shaded blue region) in a BH area Ah vs. an ADM 
angular momentum J diagram. Kerr BHs exist above the black solid line; the blue dotted line corresponds 
to the fundamental existence line for scalar clouds and the red solid line to the BS limit. We display curves 
with constant ADM mass M connecting Kerr BHs (solid green) and KBHsSH (dashed black). In the region 
of non-uniqueness - above the black solid line and below the blue dotted line - the KBHsSH always have 
larger area for the same M and J. (Right panel) Domain of existence of KBHsSH (shaded blue region) in 
a BH area Ah vs. BH temperature Th diagram. The solid green line corresponds to the Kerr limit, being 
delimited by extremal BHs (Th = 0) and arbitrarily small BHs (Th oo). Along constant Q# lines (or 
equivalently, constant w lines), KBHsSH can either get to extremal - as for D#/M = 0-9 or become arbitrarily 
small and hence high temperature at both limits of the line, cf. Figure [4] (left panel). 




Figure 6: (Left panel) Constant Q KBHsSH in a ADM mass vs. scalar field frequency diagram. (Right 
panel) Constant q KBHsSH in a J/M 2 vs. ADM mass diagram. 


4.3.2 A sample of reference configurations 

To illustrate the solutions of KBHsSH, as well as the limiting cases, we present five examples of qualitatively 
different cases, all with m = 1, n — 0, namely: 

I : a typical BS, with frequency w = 0.85; 

II : a Kerr BH in the region of non-uniqueness. We have chosen its mass and angular momentum to be 
M ~ 0.415, J ~ 0.172; this corresponds to rn — 0.066; 
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III : a KBHSH in the region of non-uniqueness with the same M, J as the Kerr BH in example II. This 
KBHSH is Kerr-like and has w = 0.975 and r# = 0.2; 

IV : a KBHSH with w = 0.82 and rn = 0.1 (this solution is close to the main branch of BSs which is the 
most stable branch); 

V : a KBHSH with w = 0.68 and r# = 0.04 (the nearby BSs are in the secondary branch, being unstable). 

In Appendix [C] we provide a set of plots that describe properties of these example solutions. Each figure 
has ten (or eight) panels divided into two colums and five (or four) rows. The figures exhibit the following 
quantities: 


• the metric functions Fo,Fi,F 2 ,IT, cf. equation (2.5), in Figures [TTj[12}[l3|[l4| respectively; 

• the metric coefficient g t t, in Figure [l5j 


the scalar field amplitude </>, cf. equation (2.6), in Figure 16 


• the scalar field energy-momentum tensor component E = — Tjf, in Figure 17 


• the scalar field energy-momentum tensor component J = X^, in Figure 

• the Ricci scalar R, in Figure [19] 

• The Kretschmann scalar K = R a bcdR abcd , in Figure 
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In all of the figures described above, Figures [IT[j20j the left column displays 3D plots, whereas the 
right column shows 2D plots of the corresponding function in terms of the radial variable for three different 
angular coordinates. Concerning the left column, the axes for the 3D plots are p = rsinO (with p > ru) and 
z = rcosO (with — < z < £ maa: , where the value of the is chosen for convenience for each case). 

In these 3D plots, only the near horizon region is shown. Concerning the 2D plots, we show the full radial 
dependence of the functions for three different angles: 0 = 0 (red line), 0 = 7r/4 (blue line) and 0 = n/2 
(green line). We recall that the solutions are invariant under a reflection on the equatorial plane. The radial 
coordinate there is x = 1 — r^/r (for BHs), and x = r/ (1 -hr) (for the BS example), such that the asymptotic 
values are approached at x = 1. 

The numerical data for these reference configurations can be found in [8l] . 


4.3.3 Horizon and ergo-regions 


As for Kerr BHs, KBHSH have a topologically spherical horizon at r = rn- Geometrically, however, the 
horizon is a squashed sphere. This can be seen by evaluating the circumference of the horizon along the 
equator, 

L e = 2tt r H e F ^ rH ^/ 2) , (4.4) 

and the circumference of the horizon along the poles, 


2 r H f d8e F2{rH ’ e) 

Jo 


(4.5) 


In Figure [7] (left panel) we show the ratio of the equatorial circumference to the polar circumference for some 
KBHsSH. As expected the squashing of the horizon produced by the rotation is such that L e /L v is typically 
larger than one; but close to the secondary branch of BSs, one finds L e /L v slightly smaller than one. 

On the horizon, the scalar field profile function 0(r#, 0), “energy density” —T t t and “angular momentum 
density” vary with the angular coordinate, as it is manifest in Figures 
solutions discussed in Section 4.3.2 In Figure [8] we show the scalar field va 
and its “angular momentum density' 


16 17 and 


18 


for the sample 


ue, its “energy density” —T\ 


T* for a KBHSH solution with r# = 0.2, m = 1 and w = 0.9. Note 
in particular that both the scalar field value and the “angular momentum density” vanish on the horizon 
poles - a property common to all solutions - but not the “energy density”. 
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Figure 7: (Left panel) The ratio L e /L v for three sets of KBHsSH solutions with fixed values of Qh (or, 
equivalently, w). In the three cases, the KBHsSH interpolate between a BS and i) another BS (£Ih/H = 0.79), 
ii) an extremal KBHSH (Q#/yU = 0.88) and iii) a Kerr BH = 0.96). (Right panel) The horizon area 

vs. q for the same three sets of solutions. 



Figure 8: The scalar field (left panel), its “energy density” (middle panel) and “angular momentum density” 
(right panel) on the horizon, for an example of a KBHSH (all panels have been multiplied by a factor of 
10 3 ). 


It would be interesting to perform a study of the Gaussian curvature of the horizon spatial sections. It is 
well known that for Kerr, beyond a certain rotation parameter, this curvature becomes negative on the polar 
caps 82 . This, in turn, prevents a global isometric embedding in Euclidean 3-space, albeit such embedding 
is possible in curved embedding 3-spaces, as hyperbolic 3-space 83 . Since an isometric embedding is a tool 


to gain insight into the geometry of BHs (see e.g. 83,84 and references therein) a study of the Gaussian 
curvature and global isometric embeddings for KBHsSH would be useful. 

We now turn our attention to the ergo-surfaces and ergo-regions of KBHsSH. These have been studied 
in detail in 52 . Three different qualitative cases can be actually seen in Figure [l5j concerning the example 
solutions of the previous subsection. For example I, g tt is always negative and there is no ergo-region. 
Indeed this is the case for some BSs. Some other BSs have a toroidal ergo-surface: an ergo-torus. Example 
II, the Kerr BH example, exhibits the usual ergo-sphere. An ergo-sphere also occurs for examples III and 
IV. Example V, exhibits a more complex structure: there is both an ergo-sphere and an ergo-torus, i.e. it 
exhibits an ergo-Saturn. In Figure [9] we exhibit the distribution of the different structure of ergo-regions in 
the M-w diagram for KBHsSH. 
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(no ergo-surface) 



Figure 9: Ergo-surface distribution in the domain of existence for KBHsSH with m = 1 and n = 0 (shaded 
blue region) in an ADM mass vs. scalar field frequency diagram. Adapted from 52 . The 3D plots illustrate 


the ergo-surfaces and, for KBHsSH, also the horizon, which is the innermost half-sphere. These are obtained 
from the numerical data for specific solutions. 


4.3.4 Extremal KBHsSH 


With the ansatz (2.5), extremal KBHsSH are found as a different th 0 limit of the hairy solutions from 
that of BSs (see Appendix [A] for an explicit proof in the Kerr limit). In the former case, however, a number 
of metric functions diverge; as a result, only near-extremal solutions can be constructed within the scheme 
described above. Therefore, extremal KBHsSH solutions are studied by constructing them directly, within 
the same numerical scheme, although for a different metric ansatz. A systematic study of these configurations 
will be presented elsewhere; here we review their basic properties only. 

Our approach to study directly extremal KBHsSH is to use a slightly different version of ( |2.5| ): again 
with four unkwnown functions (i^, W ) and the same scalar field ansatz, but now we take N = (1 — r^/r) 2 , 
where r# > 0 is a constant. Note that the extremal Kerr solution can be written in this form. The 
numerical treatment of the problem is similar to the non-extremal case; in particular, a new radial coordinate 
x = yjr 2 — r 2 H is again introduced. The boundary conditions are similar to those used in the non-extremal 
case, except for the scalar field 0, which now vanishes on the horizon. 

The extremal KBHsSH have finite horizon size and global charges and possess a regular horizon in terms 
of the Kretschmann invariant; they may, however, exhibit other more subtle pathologies, see e.g. 85 . Fixing 


(m, n), these solutions form a line in the pararameter space of KBHsSH, delimiting a part of the boundary 
of the domain of existence, as described in Section |4.3.1| This line appears to have the same qualitative 


behaviour of the BS line; however, instead of the zero mass, flat space limit, the extremal KBHsSH line 
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starts at a point corresponding to a scalar cloud line around an extremal Kerr background. Incidentally, 
we remark that the analytical estimates in 47 provide a very good approximation for the position of this 
point. Then, in a mass-frequency diagram, this line increases up to a point corresponding to a maximal value 
of mass; then it decreases until a minimal value of the frequency is approached, where is backbends and 
keeps decreasing, until a minimal value of the mass is reached. We expect this curve to inspiral towards a 
central value where, we conjecture, it meets the endpoint of the BS spiral in a singular solution for a critical 
configuration. This is the behaviour noticed for five dimensional extremal BHs with scalar hair in Anti-de 
Sitter spacetime 85 , and also (up to some details) in a Minkowski spacetime background 


5 Phenomenological properties 

If fundamental, sufficiently stable, scalar fields exist in Nature, of the type that can source KBHsSH, the 
solutions exhibited herein could play a role in astrophysical systems. Moreover, if such fields only interact 
gravitationally, as some dark matter candidates, strong gravity systems may be the only arenas where they 
leave observational signatures. As such it is of interest to understand physical properties of KBHsSH with 
phenomenological relevance. This discussion becomes even more interesting in view of the observation in [ 1 
that some of these properties deviate considerably from those of the Kerr solution, an uncommon feature 
for BHs in alternative theories of gravity (see e.g. the discussion in 87). In the following we shall expand 
on two properties of phenomenological interest, already briefly discussed in I]: the quadrupole moment 
and the orbital frequency at the innermost stable circular orbit (ISCO), and mention some other possible 
phenomenological directions of research. 


5.1 Quadrupole Moment 


In Newtonian gravitational systems, the multipolar expansion provides a complete description of the gravi¬ 
tational field for a distribution of (static) masses. In relativistic gravity, Geroch 88 and Hansen 89 have 


devised strategies to define a physical significant multipolar expansion. In particular, the quadrupole mo¬ 
ment is of great interest; because in principle it can be measured, say, by gravitational wave signals of a star 
or a smaller BH orbiting a central (larger) BH (see e.g. |90l) and it can be used to test the no-hair idea 13 


since for the Kerr solution the quadrupole moment is completely determined in terms of the ADM mass M 
and angular momentum J as —J 2 /M. 

In computing the quadrupole moment of KBHsSH, we shall follow the general procedure described 
in [90 92 for extracting it from the asymptotics of a stationary and axially symmetric spacetime. The above 
references use a line element written in quasi-istropic coordinates, with 


ds 2 = e 2 ^~ v) ( dR 2 + R 2 d6 2 ) + R 2 sin 2 0B 2 e~ 2,y (dy - Wdt) 2 - e 2lJ dt 2 


(5.1) 


where £, z/, B and W are functions of R,Q. The quadrupole moment of the spacetime is encoded in the 
large -R asymptotics of the metric functions z/, B\ to leading order, one finds 
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(with P n the Legendre polynomials) such that the quadrupole moment Q is given by 

B 0 


Q = -V2 
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(5.2) 


(5.3) 




In evaluating (5.3) for KBHsSH solutions, we need first to bring the line element (2.5) into the form 
This is done by using the coordinate transformation r = R( 1 + Jr) 2 - Then the coefficients Bo and which 
enter the asymptotics ( |5.2| ) can be read from the far field form ( |B.3 ), (B.4) of the solution (expressed in 
terms of the quasi-isotropic R ), 

B 0 = a 5 - j^( 2c t - rH) 2 , V 2 = ^ (-24&1 - 8 a 5 c t + 2c? + 4a 5 r# - 3 c 2 r H - 3c t r 2 H ) . (5.4) 
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In Figure 10 (left panel) we exhibit the reduced quadrupole (i.e the ratio of the KBHSH quadrupole to 
the Kerr quadrupole) for a set of q =constant as well as Qh =constant KBHsSH (and BSs) in terms of the 
dimensionless parameter J/M 2 . One observes that KBHsSH can have a quadrupole moment one order of 
magnitude larger than Kerr BHs within the Kerr bound and up to two orders larger than Kerr solutions 
beyond the Kerr bound. 



Figure 10: (Left panel) reduced quadrupole. (Right panel) angular frequency at the ISCO. From 


0 - 


5.2 Orbital frequency at the ISCO 

The study of geodesics represent a standard way to analyze a given spacetime geometry. Moreover, in 
astrophysical environments it is thought that the edge of accretion disks is determined by the ISC Os around 
a given BH. The accelerated charges that orbit the BH will emit synchroton radiation which, in the simplest 
model, will have a cut-off at the orbital frequency of geodesics at the ISCO. Thus, measurements of the ISCO 
via observations of accretion disks can, in principle, be used to evaluate the properties of an astrophysical 
BH. As such we have studied the angular frequency at the ISCO for a large set of KBHsSH. 

The geodesic motion is studied along the equatorial plane, 0 = tt/2; then the Lagrangian of a timelike 
test particle (the only case studied here) is 

2 C = e 2Fl T — + e 2 F 2 r 2 (ip - Wi ) 2 - e 2 F °Ni 2 = -1 . (5.5) 

Note that i^, W depend only on r; also a dot denotes a derivative w.r.t. the proper time. The stationarity 
and axisymmetry of the KBHsSH metric implies the existence of the first integrals 

e 2F2 r 2 (ip-Wt) = L , (5.6) 

(e 2F °N - e 2 F 2 r 2 W 2 )i + e 2 F 2 r 2 W<p = E , 


with E and L the specific energy and angular momentum of the test particle. Then the orbital angular 
velocity is expressed as 


n _ <P _ w e 2F °~ 2F2 LN 
c i r 2 (LW — E) 


(5.7) 


The equation governing the variation of the radial coordinate r for an orbit on the equatorial plane is 


r 2 = V(r) = e 


-2Fi 


N 


- 2 f 2 L 2 , e-^jE-LWf 


72" + 
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(5.8) 
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To determine circular orbits we need the first derivative of V, which is 


V'(r) = e~ 2Fl \2NF[ - 2 e~ 2F °(E - LW) 2 {Fq + F[) - N' 

e 2F -L 2 
d-5— 


(5.9) 


(2N(1 + r(F[ + F?)) - 2N') + 2 e~ 2F °L{LW - E)Wj . 


The second derivative of V is also of interest; however, its expression is very long and not particularly 
enlightening; so we shall not exhibit it here. 

The requirement for a circular orbit at r = r c is V(r c ) = V'(r c ) = 0. From (5.8), (5.9), this results in two 
algebraic equations for E , L which are solved analytically, possessing two distinct pairs of solutions (E + , L+) 
and (i£_,L_), corresponding to co-rotating and counter-rotating trajectories. 

The solutions for E,L are then replaced in the expression of V"{r c ), asking for V"(r c ) < 0. For the 
configurations studied so far we have noticed a (qualitative) analogy with the Kerr BH. First, the circular 
geodesic motion is possible for r c > r m i n only, a constraint imposed by the requirement for the energy E 
to be real. Then for r m * n < r c < risco only unstable circular orbit can exist, i.e. with V"(r c ) < 0. For 
r c > risco the circular orbits are stable. 

In Figure [lO] (right panel) we exhibit the angular frequency at the ISCO for co-rotating and counter¬ 
rotating geodesics and one can observe the deviations relative to the Kerr solution 11]. A similar study, but 
describing the angular frequency at the ISCO in terms of the horizon mass and angular momentum will be 
presented elsewhere 


79 


5.3 Other phenomenological studies 

As discussed in the Introduction, the next decade may bring the first detection of gravitational waves, and 
with it a tool to test the general relativity BH paradigm. It is therefore important to understand how alter¬ 
native models to this paradigm, for instance models with scalar fields, can have different phenomenological 
consequences. 

The strongest gravitational wave signal is expected to occur in the merger of two compact objects (such as 
neutron stars or BHs). To study such mergers, the only available tool is to resort to fully non-linear numerical 
relativity simulations. For the case of scalar-tensor theories of gravity, numerical relativity simulations are 
still in their infancy, e.g. 193-95], but qualitatively different features have already been reported - see 96 for 
a review. For the case of minimally coupled scalar fields, numerical simulations, and in particular collisions, 


of BSs have been performed (see 59 for a review), but not yet of KBHsSH. This will certainly be of interest, 
also to test the stability of these solutions - see Section [6j 

The interaction of minimally coupled scalar fields with BHs has, nevertheless, been considered using 
various approximation schemes, with the goal of extracting gravitational wave emission, e.g 
particular, in 


100 


97101 


In 


it was shown that the gravitational wave response of a perturbed BH surrounded by a 
‘dirty’ environment composed of a quasi-stationary scalar cloud has an imprint in its late time tail. Even 
though such tail is unlikely to be measured in the near future, this example suggests how the gravitational 
wave signal from a KBHSH can be different from that of the conventional BHs in general relativity. 

Another potentially observable feature of BHs also discussed in the Introduction is the BH shadow. 
Firstly, the presence of a shadow is regarded as a smoking gun for the presence of an event horizon (more 
precisely, of an apparent horizon). By contrast, for BSs there are unbound geodesics that can approach 
arbitrarily close to the centre of the BS 63 and there is no shadow (see also [102 ). As KBHsSH interpolate 


continuously between Kerr BHs and BSs, their shadows should vary continuously between the Kerr shadow 
and no shadow at all. It is of interest to understand quantitatively how this variation occurs and how 
distinctive can these shadows be. 
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6 Stability? 


6.1 General comments on stability of BHs 

The stability properties of any exact (analytical or numerical) solution of the Einstein equations are central 
to its physical relevance. It is important to remark, however, that absolute stability is not mandatory for 
physical relevance. Metastable states (e.g. Uranium or Plutonium) can play a physical role and, ultimately, 
essentially all systems in the Universe are only metastable. Thus, if instabilities are present, the discussion 
should focus on how large/small the timescales of the instabilities are, as compared to other typical timescales 
in the physical processes where the solution may play a role. 

Establishing the stability of BH solutions is highly non-trivial. One may address, by increasing order 
of complexity, the i) mode stability, ii) linear stability or iii) fully non-linear stability. Mode stability, in 
vacuum , was established for Schwarzschild BHs in the classical works of Regge-Wheeler 103 and Zerilli |104| , 


and for Kerr in the work of Whiting 105 . It relies on an analysis of quasi-normal modes of the BHs. In 


the Kerr case, the analysis is made possible by the algebraic specialty of the solution (Petrov type D), 
which allows the decoupling of gravitational perturbations, using the Newman-Penrose formalism. We have 
verified that KBHsSH are (generically) algebraically general (Petrov type I). Thus, no decoupling of the 
coupled gravitational-scalar system is expected to occur and thus a similar mode analysis to that performed 
for Kerr cannot be carried through. 

Mode stability does not, however, guarantee linear stability (quasi-normal modes do not form a basis 
of the space of functions). As an illustration of this difference, recently, it has been argued that extremal 
Kerr BHs are unstable against a special class of linear perturbations 1106 . Even for non-extremal Kerr, 
the linear stability has not yet been established (see e.g. [107]). Finally, linear stability does not guarantee 
non-linear stability. An example which has arisen great interest over the last few years is Anti-de Sitter space 
which has been shown, by numerical simulations, that it is non-linearly stable against BH formation 108 


Establishing non-linear stability is a mathematical tour de force ; for Minkowski spacetime this was proven by 
Christodoulou 109). No rigorous statement can be made for BHs, but the many fully-non linear numerical 


simulations performed with Kerr and Schwarzschild BHs have piled up evidence that, at least in the time 
and spatial scales probed by the simulations, no non-linear instabilities exist 96 . 


6.2 Stability of BSs 

Spinning BSs form a part of the boundary of the domain of existence of KBHsSH, cf. Figure [4j The stability 
of these objects is even more complicated to establish than for the (vacuum) Kerr BH case, even at the level 
of mode stability, as we shall discuss below. But let us start by overviewing non-spinning BSs (m = 0), for 
which, as for the case of BHs, stability has been analysed from various viewpoints: by considering either i) 
mode stability for linear perturbations, ii) fully numerical non-linear evolutions or iii) catastrophe theory. 

Historically, the stability of spherically symmetric BSs was first addressed in the late 80s, by means of 
a mode analysis in linear perturbation theory 110 112 . This approach leads to an eingenvalue problem, 
which is of Sturm-Liouville type and which determines the normal modes of the radial oscillations and their 
eigenvalues a 2 . The sign of the lowest eingenvalue cfq is crucial; if erg > 0 then the BS is stable, otherwise 
it is unstable. As already noticed above, the spherically symmetric BSs can be parametrized in terms of the 
central value of the scalar field 0o- Then one can show that the transition from stability to instability always 
occurs at critical points of the ADM mass M (or Noether charge, Q) against the value of the scalar held at 
the origin, where dM(0o)/^0o = dQ(4>o)/ = 0- The pattern of 0o i n a non-spinning BS is similar to that 
of 0 m ax seen in the inset of Figure [2] Observing also Figure [4] or Figure [l] one observes that critical points 
occur at the cusps of the BS line in the M — J diagram (Figure [4] or Figure [l]- right panel). A physical 
understanding of this fact is related to the binding energy of the BSs, E^ = /iQ — M: it becomes negative 
after the first critical point, indicating instability. 

Numerical non-linear evolutions of BSs started relatively recently - in the last 15 years - and allow the 
study of possible end-states of some initial data. These include dispersion to infinity of the scalar held, 
transition from an unstable to a stable configuration, or collapse to a BH. A recent review of these aspects is 
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given in Ref. 59 , together with a large set of reference. Here we only mention that the conclusions reached 


in this way on the stability question are in agreement with those within the linear perturbation theory. Note 
that, as shown by the ’no hair’ theorem in 80 all spherically symmetric BHs which occur as endstates of 


some initial data are described by the Schwarzschild vacuum metric. 
Catastrophe theory 


e.g. 


113-115 


provides yet another approach to investigate the stability of BSs (see 
116, 117| ). In this approach, an appropriate set of behavior variable(s) and control parameter(s) is 
chosen, and the series of solutions is constructed in terms of these. Under certain conditions, such a series 
generates a curve smooth everywhere except for certain points. According to catastrophe theory, passing 
through one of these points means changing the stability of the BS configurations. Between these points, 
the solutions form branches, sharing the same stability properties. The results obtained in this way confirm 
the previous conclusions. 

Concerning the stability of spinning BSs considerably less is known. Firstly, such solutions have not yet 
been analysed even within linear perturbation theory. This study appears challenging, since no decoupling 
of the coupled gravitational-scalar system is expected to occur. Additionally, no numerical evolutions of 
rotating BS initial data has yet been performed 59 , although simulations of BS binaries 


118 have found 


rotating BSs as a result of merger. To our knowledge, the only results on the stability of rotating BSs are 
those in 119 obtained by using catastrophe theory. As shown therein, the rotating BSs share a similar 
stability picture as the non-rotating solutions. In particular, the branch of rotating solutions between the 
vacuum M = J = 0 and the point where the global charges approach the absolute maximum is predicted to 
be stable. A more complicated picture occurs in the inspiraling region (see Figure [l]), with both stable and 
unstable branches. The BSs in that region, however, possess an ergo-region, cf. Figure [9] The existence of 
this region is expected to originate superradiant instabilities, as we address in the next subsection. 


6.3 Ergo-regions and superradiant instabilities 

Kerr BHs in the presence of massive bosonic fields are afflicted by superradiant instabilities that amplify 


sufficiently low frequency modes, extracting rotational energy and angular momentum from the BH 46,73- 
77,98 120-123 6 This instability is often regarded as a field analogue of the particle Penrose process [135 , 


which occurs due to the existence of an ergo-region. For a scalar field of mass g in the background of a Kerr 
BH with ADM mass M and angular momentum J, the time scale for the fastest growing unstable mode 
is 52,75 , for gM <C 1, 

1 M 2 * M 

r = —— 0.12 ms . (6.1) 

(gM) 9 J M© v ' 

As an estimate, for Mg ^0.1 and an extremal Kerr BH, J = M 2 (to minimize the time scale of the 
instability), one gets r ^ 1 day to r ~ 10 6 years, for a BHs of Mm 1M 0 to M — 1O 9 M 0 , respectively. As a 
comparison - for astrophysical relevance - the latter timescale is approximately that of the lifetime of blue 

is that for gM 


supergiant stars. An estimate from the results of 77 
two orders of magnitude. 


1 these timescales may decrease by 


KBHsSH always have an ergo-region 52 and thus should be afflicted by superradiant instabilities in the 
presence of massive bosonic fields. The same applies to BS that have an ergo-torus 65,136 . The novelty, 


with respect to Kerr, is that a massive scalar field is already present in the background solution, and thus 
one needs not consider extra fields in order to obtain superradiant instabilities, as it must be done for Kerr. 
Thus, the superradiant instabilities of KBHsSH should be seen within linear perturbations of the solutions in 


the theory (2.1). As we have argued in Section 6.1 the relevant question is the time scale of these instabilities, 


which, at the moment is not known, and a linear perturbation analysis faces the same challenges as for the 
case of BSs, discussed in Section |6.2| Let us therefore remark on some expectations, mainly comparing to 
the timescale of the instabilities observed above for the Kerr case: 


6 Superradiant instabilities should be d istinguished from superradiant scattering , which can occur for bosonic massive or 

massless fields around Kerr BHs 124 125]. Superradiant scattering of charged boso nic fie lds can also occur around spherically 

symmetric charged BHs, leading to the extraction of Coulomb energy and charge 126]. But superradiant in stabilities only 
occur around spherically symmetric charged BHs if these are enclosed in a cavity, in Anti-de-Sitter spacetime 127- 133] or in 

higher dimensional brane-world scenarios 1134 . 
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1) Superradiant instabilities of KBHsSH only occur for sufficiently high m modes; these are expected to 
be longer-lived than the fastest growing mode in the Kerr case. Indeed, as pointed out in |l| (see 
also 185]) a KBHSH solution with a given m azimuthal harmonic index should only be unstable against 
perturbations with higher azimuthal harmonic index. In the Kerr case, it is known that the timescale 
of the instabilities decays with increasing angular harmonic index £ of the spheroidal harmonic per¬ 
turbation Sim. Since a higher m requires a higher £ we expect that the unstable modes for KBHsSH 
have a larger time scale than the fastest growing mode {£ — 1) in the Kerr case. 


2) Ergo-regions of Kerr BHs are ‘larger’ than those of KBHsSH for comparable solutions, which suggests 
longer timescales for the instabilities in the latter case. In 52 it was observed that a measure of the 


size of the ergo-region (therein dubbed ergo-size ) can be introduced, which is positively correlated to 
the strength of superradiant instabilities, at least in the regime of validity of formula ( |6.1| ). Then, 
comparing the ergo-size of Kerr BHs and KBHsSH in the region of non-uniqueness, for the same global 
charges, it was observed that the ergo-size is smaller for the latter, thus suggesting instabilities are 
weaker. The region where this reasoning applies is close to BSs that have no ergo-region - see Figure [9] 
If one regards KBHsSH as bound states of Kerr BHs with BSs, it is therefore natural that these bounds 
states have a smaller ergo-region. 


3) 


In the q 
argued in 


0 limit, scalar clouds should be dynamical attractors, rather than unstable solutions 
50 that scalar clouds around Kerr BHs (i.e. the q 


It was 

0 limit of KBHsSH) are dynamically 


stable configurations; they are in synchronous rotation with the BH in a tidal locking configuration with 
some analogy to the Earth-Moon orbital tidal locking. Slightly superradiant modes are under-spinning 
whereas slightly decaying modes are over-spinning and, in both cases, dynamics tries to synchronize 
the horizon angular velocity with that of the scalar perturbation]^] This motivates from a different 
direction why a superradiant perturbation in the background of a KBHSH grows more slowly than in 
the background of Kerr: as it grows, the perturbation extracts energy and angular momentum from 
the background BH which decreases its angular velocity. Then, the background scalar mode becomes 
over-spinning compared to the horizon and competes with the under-spinning perturbation, slowing 
down (stopping?) its growth. 


To summarize, whereas spinning BSs seem to possess a branch of stable solutions, all KBHsSH should 
be afflicted by instabilities of superradiant type. It is reasonable to expect, however, that the timescales 
for (at least some) KBHsSH can be larger than those observed for Kerr solutions, which, in some cases can 
already be interesting astrophysically. Some light into this question may be shed by numerical simulations, 
for which the technology exists for evolving spinning BSs and KBHsSH 96 . We anticipate that, at least for 
a subset of KBHsSH solutions, these will appear stable even in long term numerical time evolutions, as the 
timescale for the potential instabilities will be too small. 

As a final note, let us recall that, in the region of non-uniqueness, the entropy of KBHsSH is larger than 
that of a Kerr solution with the same global charges. This prevents the former from decaying into the latter, 
adiabatically. 


7 Conclusions and Generalizations 


The main purpose of this work was to provide a detailed description of the construction and of the (basic) 
physical properties of a new type of hairy BH reported in 11 . Our results show that stationary scalar field 
matter distributions surrounding rotating BHs can exist even in the absence of scalar self-interactions or 
non-minimal couplings between the scalar field and the geometry. 

The existence of KBHsSH provides an unexpected connection between the scalar superradiant instability 
of Kerr BHs and the BSs of soliton physics. Moreover, they clarify a long standing puzzle: all field theory 
horizonless, particle-like solutions were known to possess BH generalizations except BSs. It is now clear that 
BSs also admit a BH generalization - but only if they are spinning. 


7 Bear in mind this reasoning is simplistic as the role of accretion is not properly considered; see 


101 


26 







There are many possible continuations of this work - a few already mentioned along the paper. The most 
pressing one is to deepen the issue of stability, for which a first doable step is to perform a fully dynamical 
evolution of KBHsSH. Equally relevant, especially in view of the experiments discussed in the Introduction, 
is to further explore astrophysical signatures: e.g. to compute shadows of KBHsSH, contrasting them with 

A related research line 


those for Kerr BHs, and analysing possible gravitational wave signatures 1100 . 
is to understand the dependence of phenomenological aspects on the horizon (Komar) mass and angular 
momentum rather than on the ADM global quantities. On a more theoretical side, discussing the intrinsic 
BH horizon geometry, the extremal limit and its properties and eventually considering the BH interior, seem 
worthwhile problems. 

In parallel to the aspects discussed in the previous paragraph, one may search for new solutions of BHs 
with scalar hair in more general theories, which belong to the same family as the KBHsSH considered here, 
in the sense of relying on the syncrhonization condition c a = mfln (or appropriate extensions). Indeed, the 
study in this paper, as well as those in Refs. p~[[5T[[52 , has been restricted to the simplest case of a non-self 
interacting scalar field. It should be possible, however, to adapt the approach described in Sections [2] and [3] 
to include a scalar field potential. Apart from the theoretical and technical relevance, such solutions could 
be interesting from yet another point of view. As noticed above, the maximal values of the mass of the 
BSs set an upper bound for the mass of KBHsSH, which is of the order of Mpjfi. Thus these KBHsSH 
can only be relevant in an astrophysical context if extremely light bosons exist (although they might affect 
the dynamics of small primordial BHs). The situation is very different for a self-interacting scalar field. As 
shown in 137 , when switching on a quartic self-interaction term for T, the maximal mass of stable BSs is 
of the order of A Mpjg 2 , with A the scalar field self-coupling. Therefore we expect the typical masses for 
the corresponding hairy BHs to be much larger than in the A = 0 case. 

A more involved picture is expected to exist for a complex massive scalar field with a non-renormalizable 
self-interaction potential allowing for finite mass solitons even in the absence of gravity - Q-balls. The 
existence of BH generalizations of the spinning gravitating Q-balls has been discussed in 138 at the probe 


level, i.e. for a fixed Kerr BH background. The corresponding solutions have been dubbed therein Q —clouds. 
Such configurations are also in synchronous rotation with the BH horizon, satisfying the condition w = mfln- 
In constrast to the non-selfinteracting case, however, Q-clouds exist on a 2-dimensional subspace of the Kerr 
parameter space, delimited by a minimal horizon angular velocity and by the corresponding m = l existence 
line, wherein the nonlinear terms become irrelevant and the Q-clouds reduce to the linear clouds discussed 
in Section 4.2 This implies that a more involved picture will be found when including backreaction on the 


spacetime geometry, with a much larger range for the masses of hairy BH as compared to the case discussed 
here. Furthermore, some basic features of the latter solutions should be preserved when replacing the Q-balls 
with vortons. These are a special class of scalar solitons made from loops of vortices, which are balanced 
against collapse by rotation, being the four dimensional field theory analogues of the higher dimensional 
black rings of vacuum general relativity 139 . Some of these solutions were shown to be stable 140 . We 


expect that all gravitating vortons in 141 to possess BH generalizations with many similar properties to 
the KBHsSH in this work. This type of study can moreover be used to set constraints on the properties of 

see e.g 
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the scalar field, by comparing astrophysical observations with the set of hairy BHs predictions 
for a similar approach in the spherically symmetric solitonic limit of solutions. 

Yet another possible generalization for future research is to include a matter content more generic than 
scalar fieds in the theory, and look for the corresponding hairy solutions. One could test the conjecture, 
put forward in 51 that “A (hairless) BH which is afflicted by the superradiant instability of a given field 
must allow hairy generalizations with that field. ^ In particular, this suggests the existence of Kerr BHs 
with Proca hair - see 142} 145] for superradiant instabilities of Proca fields around BHs. The simplest 
case in this line of extensions, however, are generalizations of the KBHsSH with a 17(1) gauged scalar field 
leading to Kerr-Newman BHs with scalar hair (KNBHsSH). The condition ( |3.9| ) is replaced in this case 
with w = uiQh + q&H, where q is the gauge coupling constant and <f># is the electrostatic potential on the 
horizon. Similarly, to the case described here a set of KNBHsSH emerges as backreacting charged scalar 
clouds around Kerr-Newman BHs |49l|50|. An even more complicated picture should exist when allowing for 


8 Of course, the test field must source a time and azimuthal independent energy-momentum tensor, which, for instance, 
immediately excludes a single real scalar field. 
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scalar multiplets gauged with respect to a non-Abelian gauge group. The simplest set of such solutions has 
been already discussed in a different context in 146 , for a triplet Higgs field and an SU (2) gauge field. 

Of a more theoretical interest, one may consider generalizations of the solutions in this work to different 
spacetime dimensions, D ^ 4, and also possibly with different asymptotics. Together with the self-interacting 
case, knowledge of these solutions may lead to valuable insights into the more relevant D = 4 asymptotically 
flat solutions, by establishing which properties of hairy BHs are generic, and also which ingredients are crucial. 
Some results in this directions exist already in the literature. For example, a family of asymptotically Anti-de 
Sitter rotating BHs with scalar hair and a regular horizon has been studied in 85 within D = 5 Einstein’s 


gravity minimally coupled to a complex scalar field doublet; their asymptotically flat counterparts have been 
studied in 86 . While the properties of the solutions in 85 are rather similar to the case in this work - and 


indeed was the first example of this type of solution - some of the features in 86 could hardly be anticipated. 
For example, the hair of the Dm 5 asymptotically flat hairy BH solutions is intrinsically nonlinear, since a 
Myers-Perry BH background does not allow for scalar clouds of massive, test scalar fields. 

We close with a natural question: can there be a solution of KBHsSH in an analytic closed form? In 
this context, we remark that in more that forty years, no closed form expression could be found even for the 
simpler case of spherically symmetric BSs - see however 147 for an exception in D = 3. One hope is to 


consider solution generating techniques, or in the context of supergravity models, backgrounds preserving 
some amount of supersymmetry. In the latter context, however, it is known that the existence of Killing 
spinors is incompatible with an ergo-region 148 -150 . Thus, even if hairy BHs/BSs could be found in this 


setup, they will hardly be representative of the physical properties of the generic case. 
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A New coordinates for Kerr 


As we described in the main text we have used the metric ansatz (2.5). The Kerr metric can also be written 


in this form. The corresponding expressions of the metric functions read: 


e 2Fl = (1 - — ) 2 + c t (c t - r H ) 


cos 2 9 


e 2F 2 = e -2 Fl 




F 0 = -F 2 , W = e-* Fl+F ^y/ct(ct - r H ){r H - 2c t )^p- 


(A.l) 


Expressed in this form, the solution contains two constant r# and c t . While r# fixes the event horizon 
radius, the second constant, c t < 0 does not have a transparent meaning; however, it can taken as a measure 
of non-staticity, since c t = 0 is the Schwarzschild metric. 

The expressions for various quantities of interest read (we recall G = 1 = c) 


M = ~(r H - 2ct), 

J = 7} Vct(ct^nijirH - 2ct), 
A h = 47r (r H - c t )(r H - 2 c t ), 
T h = TH 

0.H = 


(A.2) 


47T (r H ~ c t )(r H ~ 2c t y 
- r H ) 


{r H ~ c t )(r H ~ 2 ct) 

Note that the formal limit r# = 2 c t is just flat space expressed in an unusual coordinate system. 


The relation between the radial coordinate r in (A.l) and i?, the radial coordinate of the Kerr metric in 
BL form is 


r = R — 


Rh 


where 


R h = M + VM 2 - a 2 


is the (outer) horizon radius for the Kerr metric in BL coordinates and a = J/M. We notice also the simple 
relation 


th = Rh ~ 


Rh 


Observe that for extremal Kerr r# —>> 0. The coordinates 6 , ip and t are the same for both parametrizations. 


B Asymptotic expansions 

An aproximate form of the solution can easily be constructed on the boundary of the domain of integration. 
Of special interest is the expansion as r oo, which is used to evaluate the expression of the quadrupole 
moment. 

Starting with the scalar field, its asymptotic behaviour is of the form 

p -\/ Li 2 ~w 2 r 

<Kr,0)=m -+ •••. (B.l) 
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as to describe bound state solutions of the KG equation, which requires w < fi. The function f(0) which 
enters the asymptotics of the scalar field can be expressed as 


f(e) = J2fkP™+2k(cose), 


(B.2) 


k=0 


where P™+ 2 k ( c °s 0) are the associated Legendre functions and /& are constants fixed by numerics. 

One can also construct an approximate expression for the metric functions in inverse powers of r (note 
that for the purposes here the scalar field contribution can be neglected, since </> decays faster than any 
power of r). The leading order terms in such asymptotic expansion read 

F lr a) _ °t , °t r H MO) 

F o(r,V) - j + ^~ + — 3 - + ..., 


f 2 (r> o) = ~* Ct{ct + Th) ■ f23{d) 

W(r,6) = ^ + ^ + ^l 


(B.3) 


where the expressions for the ^-dependent coefficients are: 

fo3(0) = h + - I (246 a + 8 a 5 ct - 2 c 3 - 4a 5 r H + 3c 2 r H + 3 c t r 2 H ) cos 2 0, 

fn(9) = ~j(ct + r H ) + a 5 cos20, 


fi3(0) = ^ ( 8&1 + “ 5 ( 8 c t - 4 r H ) - c t {2c t + c t r H + r H ) 

+ (246i + 8a^c t — 2cj* + Yla^ru + 3r# + 3 c t r 2 H ) cos 2# 

/23(0) = ^7 ^ 86 i + 4a 5 (2c t + 3 r H ) - c t (2c 2 + c t r H + r 2 H ) (B.4) 

+ (24&i + 8 a 5 c t — 2(? t — ^cl^th + 3 c^vh + 3 c t r 2 H ) cos 2# 

n; 5 ( 6 >) = ^ ^ - 36a 5 c^ + 15c^ct(7ct + r H ) + 45n; t + 75n;t cos 20^j . 

In these expressions, q, as, fq and w t are arbitrary constants which are extracted from the numerical 
output. Note that for the Kerr BH, the solution have 


w t 


2 _ \ 

= —c t v / c t (c t - r H )(ct - r ff )(2ct - r^), 61 = -c 2 (2c t - 3r#), 


15 

1 


«5 = jct(c t - r H ), c v = sfcjct - - 2 c t ). 


C Reference solutions plots 


In this appendix we provide the plots mentioned in Section 4.3.2 
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Figure 11: The metric function Fq for example solutions I-V, cf. Section 4.3.2 
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Figure 12: The metric function F\ for example solutions I-V, cf. Section 4.3.2 
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Figure 13: The metric function F 2 for example solutions I-V, cf. Section 4.3.2 
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Figure 14: The metric function W for example solutions I-V, cf. Section 4.3.2 
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Figure 15: The metric coefficient g tt for example solutions I-V, cf. Section 4.3.2 
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Figure 16: The scalar field profile function (j) for example solutions I and III-V, cf. Section 4.3.2 
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Figure 17: The scalar “energy density” for example solutions I and III-V, cf. Section 4.3.2 
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Figure 18: The scalar “angular momentum density” for example solutions I and III-V, cf. Section 4.3.2 
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Figure 19: The Ricci scalar for example solutions I and III-V, cf. Section 4.3.2 
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Figure 20: The Kretschmann scalar for example solutions I-V, cf. Section 4.3.2 
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